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Scaling behavior of quark propagator in full QCD 
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We study the scaling behavior of the quark propagator on two lattices with similar physical volume 
in Landau gauge with 2+1 flavors of dynamical quarks in order to test whether we are close to the 
continuum limit for these lattices. We use configurations generated with an improved staggered 
("Asqtad") action by the MILC collaboration. The calculations are performed on 28 3 x 96 lattices 
with lattice spacing a — 0.09 fm and on 20 3 x 64 lattices with lattice spacing a — 0.12 fm. We 
calculate the quark mass function, M(q 2 ), and the wave-function renormalization function, Z(q 2 ), 
for a variety of bare quark masses. Comparing the behavior of these functions on the two sets of 
lattices we find that both Z(q 2 ) and M(q 2 ) show little sensitivity to the ultraviolet cutoff. 

PACS numbers: 12.38.Gc 11.15.Ha 12.38.Aw 14.65.-q 



I. INTRODUCTION 

Quantum Chromodynamics (QCD) is widely consid- 
ered to be the correct theory of the strong interactions. 
Quarks and gluons are the fundamental degrees of free- 
dom of this theory. The quark propagator contains valu- 
able information about nonperturbative QCD. The sys- 
tematic study of the quark propagator on the lattice 
has provided fruitful interaction with other approaches 
to hadron physics, such as instanton phenomenology 0, 
chiral cjuark models 0] and Dyson-Schwinger equation 
studies 0,3. As a first principles approach lattice QCD 
has provided valuable constraints for model builders. In 
turn, such alternative methods can provide feedback on 
regions that are difficult to access directly on the lattice, 
such as the deep infrared and chiral limits. 

The quark propagator haspreviously been studied us- 
ing Clover 

Ull, staggered 0J1 and Overlap UEEl 



actions. For a review, see Ref. |12(. All these actions have 
different systematic errors and the combination of these 
studies has given us an excellent handle on the possible 
lattice artifacts in quenched QCD. 

In this study we focus on the Landau gauge quark 
propagator in full QCD, and extend our previous 
work 1 1311 to a finer lattice with lattice spacing a = 0.09 
fm |14| but similar physical volume in order to test 
whether we are close to the continuum limit for these 
lattices. The scaling behavior of the momentum space 
quark propagator is examined by comparing the results 
on these two lattices. Our results show that there are no 
significant differences in the wave-function renormaliza- 
tion function and quark mass function on the two sets of 
lattices. Therefore the scaling behavior is good already 
at the coarser lattice spacing of a = 0.12 fm. 

The configurations we use in this study were gener- 
ated by the MILC collaboration 114 ^. l3~5i and are avail- 
able from the Gauge Connection jig. The dynamical 
configurations have two degenerate light fermions for the 



u and d quarks and a heavier one for the strange quark. 
Weighting for the fermion determinants is provided by 
the so-called, "fourth root trick." . While the current nu- 
merical results 17] provide compelling evidence that the 
fourth root trick gives an accurate estimate of the dy- 
namical fermion weight, the formal issue of proving that 
this provides the determinant of a local fermion action 
from first principles remains unresolved. 



II. DETAILS OF THE CALCULATION 

The quark propagator is gauge dependent and we work 
in the Landau gauge for ease of comparison with other 
studies. Landau gauge is a smooth gauge that preserves 
the Lorentz invariance of the theory, so it is a popu- 
lar choice. As derived in Ref. an improved Landau- 
gauge-fixing functional, 
where 
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2 is used 



F?[{U}\ = £ ~Tr {U<j(x) + U^( X y} , (1) 
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^ = E \^{U°{x)UC(x + ,i) + h.c.} 



X , (1 



U°(x) = G(x)U lt (x)G(x + [i)\ 



G(x) 



exp j-i£V(:r)T a j , 
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and Mo is the plaquette measure of the mean link. We 
adopt a "steepest decents" approach. The functional 
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derivative of with respect to u a provide 

M*) = E l u »(? - /*) - M*) - h.c] tracclcss (5) 



C/ M (a; - 2ii)U ll (x - fx) 



-UJx)UJx + ju) - h.c. 



(6) 



tracclcss 



and 



A Imp (z) = -Ai(i) - -^A 2 (z). 
The resulting gauge transformation is 

Gimp (as) = exp|^Ai mp (x)| , 



(7) 



(8) 



where a is a tuneable step-size parameter. The gauge 
fixing algorithm proceeds by calculating the relevant Aj 
in terms of the mean-field-improved links, and then ap- 
plying the associated gauge transformation, Eq. (JSJ, to 
the gauge field. The algorithm using conjugate gradient 
Fourier acceleration is implemented in parallel, updating 
all links simultaneously, and is iterated until the Lattice 
Landau gauge condition 



^Imp 
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Tr {Ai mp (a;)Ai m p(a;) t } 



(9) 



is satisfied with accuracy of 6i < 10~ 12 . 

As this gauge fixing finds a local minimum of the gauge 
fixing functional, we are necessarily sampling from the 
first Gribov region. Our ensemble contains no gauge- 
equivalent configurations and hence has no Gribov copies 
as such. However, our configurations are local minima 
and absolute minima and therefore are not from the Fun- 
damental Modular Region 19]. It is known from pre- 
vious SU(3) studies that neither the gluon nor quark 
propagator display any obvious Gribov noise above and 
beyond the ensemble statistical noise and so we do not 
consider it further here It will be interest- 

ing to repeat this calculation for the Gribov-copy free 
Laplacian gauge, and to do a systematic search for Gri- 
bov noise in Landau gauge , but these are left for future 
studies. 

The MILC configurations were generated with 
the 0(a 2 ) one- loop Symanzik- improved Liischer-Weisz 
gauge action m. Th e dy namical configurations use the 
Asqtad quark action |2J], an 0(a 2 ) Symanzik- improved 
staggered fermion action which removes lattice artifacts 
up to order a 2 g 2 . We refer to the a = 0.09 fm lattice as 
the "fine" lattice and the a — 0.12 fm one as the "coarse" 
lattice. 

We explore two light sea quark masses, ma — 0.0062 
(to = 14.0 MeV) and ma = 0.0124 (to = 27.1 MeV). The 



TABLE I: Lattice parameters used in this study. The dy- 
namical configurations each have two degenerate light quarks 
(up/down) and a heavier quark (strange). The light bare 
quark masses for the 28 3 x 96 lattice are 14.0 MeV and 27.1 
MeV with a strange quark mass of 67.8 MeV. For the 20 3 x 64 
lattice the bare quark masses range from 15.7 MeV to 78.9 
MeV. The lattice spacing is a ~ 0.12 fm for the 20 3 x 64 lattice 
and a ~ 0.09 fm Q for the 28 3 x 96 lattice. 





Dimensions 


j3 a Bare Quark Mass 


# Coring 


1 


28 3 


x 96 


7.09 0.086 fm 14.0 MeV, 67.8 MeV 


108 


2 


28 3 


x 96 


7.11 0.086 fm 27.1 MeV, 67.8 MeV 


110 


3 


20 3 


x 64 


6.76 0.121 fm 15.7 MeV, 78.9 MeV 


203 


4 


20 3 


x 64 


6.79 0.121 fm 31.5 MeV, 78.9 MeV 


249 


5 


20 3 


x 64 


6.81 0.120 fm 47.3 MeV, 78.9 MeV 


268 


6 


20 3 


x 64 


6.83 0.119 fm 63.1 MeV, 78.9 MeV 


318 



bare strange quark mass was fixed at ma — 0.031, or to = 
67.8 MeV for a — 0.09 fm. The values of the coupling 
and the bare light sea-quark masses are matched such 
that the lattice spacing is held constant. The simulation 
parameters are summarized in Table [I] with the lattice 
spacings taken from [T3 ]. 

On the lattice, the bare propagator S(a;p 2 ) is related 
to the renormalized propagator S ren (/j,;p 2 ) through the 
renormalization constant 



S(a;p 2 ) = Z 2 (a;riS ren (^p 2 ). 



(10) 



In the continuum limit, Lorentz invariance allows one to 
decompose the full quark propagator into Dirac vector 
and scalar pieces 



S- 1 (p 2 ) = Z- 1 (p 2 )[ il -p + M(p 2 )}, 



(11) 



where M(p 2 ) and Z{p 2 ) are the nonperturbative mass 
and wave function renormalization functions, respec- 
tively. Asymptotic freedom implies that, as p 2 — ► oo, 
S(p 2 ) reduces to the tree-level propagator 

S^V) -^ij-p + m, (12) 

up to logarithmic corrections. The mass function M is 
renormalization point independent and for Z we choose 
throughout this work the renormalization point as /i = 
3.0 GeV, i.e., 



,,,,, . S(a;/i 2 ) 



Z 2 (a;n) 



= 1, 



(13) 



thus defining Zi(a\ /x). 

The tree-level quark propagator with the Asqtad ac- 
tion has the form 



S X (p) = «E 7 m9(Pm) 



(14) 
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FIG. 1: The unquenched wave-function renormalization func- 
tion Z(q 2 ) and mass function M(q 2 ) for a variety of valence 
quark masses (shown in the inset), with the light sea-quark 
mass fixed at m = 14.0 MeV. The renormalization function 
is renormalized at q = 3.0 GeV. 



where q(Pn) is the kinematic momentum given in 7] 

g M = sin(p M )[l + - sin 2 (p M )] . (15) 

The 7 M form a staggered Dirac algebra (see Eq.(A.6) and 
(A.7) of Ref. 0). Having identified the kinematic mo- 
mentum, we define the mass and renormalization func- 
tions by 

S-'ip) = Z-\q) [»£ (7» M ) + M(q)] . (16) 
Additional details can be found in Ref. 0] 

III. NUMERICAL RESULTS 

In Fig. n we show the results for the mass function 
M(q 2 ) and wave- function renormalization function Z(q 2 ) 
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FIG. 2: The unquenched wave-function renormalisation func- 
tion for the two different values of the light sea quark mass 
on the fine lattice (14.0 MeV and 27.1 MeV). The valence 
quark masses are m — 14.0 MeV (top) and m = 135.6 MeV 
(bottom), the lightest and heaviest in our current sample re- 
spectively. The renormalization function is renormalized at q 
= 3.0 GeV. 



for the lightest of our light sea quark masses for a variety 
of valence quark masses. In these figures, one valence 
quark mass (14.0 MeV) is identical to the light sea quark 
mass, as in full QCD. The others are partially quenched 
results. The data are ordered as we expect, i.e., the 
larger the bare valence quark mass, the higher M(q 2 ). 
The wave- function renormalization function, Z{q 2 ), on 
the other hand, is infrared suppressed and the smaller 
the valence quark mass, the more pronounced the dip 
at low momenta. In Figs. |21 and |3] we instead hold the 
valence quark mass fixed and vary the sea quark mass. 
Clearly the dependence over this small range of sea quark 
masses is weak. Unfortunately we only have two dynam- 
ical sets to compare, and for the lightest valence quark 
the data are rather noisy. 

Next we work on two lattices with different lattice 
spacing but similar physical volume. We compare the 
wave- function renormalization function Z(q 2 ) and mass 
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function M(q 2 ) for two lattices with different lattice spac- 
ing a in full lattice QCD. 

In Fig. we show the quark propagator from the fine 
lattice for full QCD (light sea-quark mass and valence 
quark mass equal) with the light quark mass set to m = 
27.1 MeV. This is compared with data from the coarse 
lattice by a simple linear interpolation from the four dif- 
ferent data sets so the running masses are the same at q 2 
= 3.0 GeV. Fig. repeats this for the lighter sea quark, 
m = 14.0 MeV. The quark propagators are in excellent 
agreement, showing no dependence on the lattice spac- 
ing. 



IV. CONCLUSIONS 

In this study we performed a systematic comparison of 
the Asqtad quark propagator in full QCD for two lattices 
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FIG. 4: Comparison of wave-function renormalization func- 
tion Z(q 2 ) and mass function M(q 2 ) for two different lat- 
tices. Triangles correspond to the quark propagator at mass 
m = 27.1 MeV from 28 3 x 96 with lattice spacing a = 0.09 
fm. The open circles are the data from 20 3 x 64 with lattice 
spacing a = 0.12 fm obtained by interpolating four different 
set of light quark masses making the M(q 2 ) value matched 
for both lattices at q = 3.0 GeV. The renormalization point 
for Z(q 2 ) is set at q = 3.0 GeV for both lattices. 



with different lattice spacing in order to establish how 
close these lattices are to the scaling region and hence 
to the contiuum limit. We compared the two functions 
Z(q 2 ) and M(q 2 ) on fine and coarse lattices and found 
them to be consistent within errors. We can thus deduce 
that for both lattices we are close to the scaling region 
for the quark propagator, which for example makes these 
lattices suitable for future studies attempting to deter- 
mine quark masses |25j| . 
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FIG. 3: The unquenched quark mass function for the two 
different values of the light sea quark mass on the fine lattice 
(14.0 MeV and 27.1 MeV). The valence quark masses are 
m = 14.0 MeV (top) and m = 135.6. MeV (bottom), the 
lightest and heaviest in our current sample respectively. 
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FIG. 5: This figure is same as figure 3, except the light quark 
mass of 28 3 x 96 with lattice spacing a = 0.09 fm is m — 14.0 
MeV. The renormalization point for Z(q 2 ) is set at q — 3.0 
GeV for both lattices. 
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